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The construction of the southern site of the Pierre Auger Observatory is almost completed. 
Three independent measurements of the flux of the cosmic rays with energies larger than 
lO'^* eV have been performed during the construction phase. The surface detector data 
collected until August 2007 have been used to establish a flux suppression at the highest 
energies with a 6cr significance. The observations of cosmic rays by the fluorescence detector 
allowed the extension of the energy spectrum to lower energies, where the efficiency of the 
surface detector is less than 100% and a change in the spectral index is expected. 



1 Introduction 

Cosmic rays are particles that travel through the galactic and intergalactic space, arriving on 
Earth in a broad energy range up to 100 EeV. Their flux drops steeply with energy, from a few 
particles per second per m^ up to one particle per km^ per century at the highest energies. 
The shape of the flux depends on the evolution of the sources, the mechanisms that accelerate 
particles up to highest energies and the energy losses during the propagation of the particles 
from the sources to earth. Th e tra nsition from galactic to extragalactic source is expected to be 
in the 0.1-5 EeV energy range. 

In the highest energy rang e, ab ove 60 EeV, a flux suppression is expected due to the Greisen- 
Zatsepin-Kuzmin effect (GZK)'^'^ and/or due to the maximum energy that a cosmic ray acceler- 
ator can reach. The GZK suppression is a propagation effect: protons interact with the cosmic 
microwave background (CMB) radiation losing about 15% of their energy at each encounter 
with these CMB photons. Heavier elements are dissociated through photo-disintegration. 

The previous experiments agasaE] and HiRes^ have given contradictory results regarding 
the ultrahigh energy end of the energy spectrum where only indirect observations of cosmic rays 
are possible by the detection of extensive air showers. The previous measurements of the comic 
rays flux were dominated by statistical or systematic uncertainties. In order to decrease the 
uncertainties the Pierre Auger Observatory was built as a hybrid detector, combining the two 
techniques employed by the forerunner experiments: a surface detector array and a fluorescence 
detectors. Due to the hybrid technique the nearly calorimetric estimation of the energy of the 
primary particle as obtained from the fluorescence technique can be transfered to the large 
number of events recorded by the surface detector. Recently the flux suppression has been seen 
both by the Pierre Auger Observatory^ and by the HiRes collaboration^. 

The observatory is described in the second section. The surface detector data are used to 
deduce the energy spectrum above 3 EeV where the trigger efficiency is 100%, as is described 
in the third section. Another complementary data set is delivered by the fluorescence detector 
itself. It can be used to extend the energy range down to 1 EeV. This measurement is described 
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Figure 1: A typical golden hybrid event reconstruction, with an energy of 30 EeV and an incoming direction of 27° 
(left hand side) Lateral distribution. Filled circles represent acquired signals, triangles are functioning stations 
without signal used with Poisson probabilities in the maximum likelihood fit. S(IOOO) is marked with a cross, 
(right hand side) Longitudinal profile: energy deposit in the atmosphere as a function of the slant depth. 



below in the last section together with the method of combining it with the surface detector flux 
estimate. 

2 The Pierre Auger Observatory 

After entering the atmosphere, cosmic rays interact with nuclei in the air and start creating 
extensive air showers. The charged particles that reach the ground are detected with the surface 
detector (SD), their lateral spread from the air shower axis at primary energies above 10^^ eV 
is of the order of a few kilometers. On the way through the atmosphere charged particles 
excite nitrogen molecules, which afterwards emit fluorescence light in the ultra-violet band. The 
amount of light is proportional to the energy deposited by the air shower in the atmosphere. 

The Pierre Auger Observatory, located in the province of Mendoza (Argentina), is utilized 
to measure the properties of extensive air showers by observing their longitudinal development 
in the atmosphere as well as their lateral spread at ground level. The Observatory consists of 
more than 1600 water-Cherenkov detectors, filled with 12 tonnes of water each and equipped 
with three photomultipliers to detect secondary photons and charged particles. The tanks are 
spread over about 3000 km^ on a triangular grid of 1.5 km spacing. The atmosphere above 
the array is viewed by four fiuorescence detectors (FD), each housing six telescopes, located on 
the border of the area. The field of view of each telescope is 30° in azimuth, and 1.5 — 30° in 
elevation. Light is focused with a spherical mirror of 11 m^ effective area on a camera of 440 
hexagonal pixels. Each pixel is a photomultiplier tube with 18 cm-^ detection area. More details 
on detector setup and calibration can be found in^^^^. An extension of the Observatory has been 
started with a denser array of tanks equipped with muon counters which will lower 

the trigger threshold energy for the SD, and HEAT^J, three telescopes that will increase the 
field of view of FD up to 60° . The counterpart of the Southern side is in the planning phase in 
the Northern hemisphere, in Lamar, Colorado, which will provide large statistics above 50 EeV. 

An example of a reconstruction of the same air shower with the SD and FD is shown in 
Fig. [TJ The signals recorded in the tanks are converted in terms of vertical equivalent muons 
(VEM). One VEM represents the average of the signals produced in the 3 PMTs by a vertical 
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Figure 2: (left hand side) S(IOOO) attenuation in the atmosphere. The line represents an empirical fit that is used 
for the conversion to Sag.^ (right hand side) The number of events as a function of S(IOOO) (upper panel) and as 
a function of Ssg (lower panel) for equal cos^(S) intervals. 



muon that passes centrally through the SD detector. The air shower axis is obtained from the 
arrival time of the first particles in each detector station. The impact point on ground and the 
lateral distribution of signals are obtained in a global maximum likelihood minimization which 
accounts for the station trigger threshold and the overflow of the FADCs counts in the stations 
very close to the shower axis. The effect of the fluctuation of the lateral distribution function 
is minimized at 1000 m. This optimal distance is influenced by the array spacing. The signal at 
this specific distance, S(IOOO), is used as energy estimator^. 

About one in ten air showers that reach the surface detector are also observed with the 
fluorescence detector (the fluorescence detector operates only in moonless clear nights). The 
longitudinal profile of the air shower, i.e. the energy deposit as a function of traversed matter 
in the atmosphere is obtained taking into account the fluorescence and Cherenkov light contri- 
butions and the light scattering and attenuation^!^. Due to the limited field of view, the entire 
longitudinal profile usually is no recorded, so a fit with a Gaisser-Hillas function is employed to 
obtain the full profile. The energy of the cosmic ray is the integral over this function with a 
correction of (1 ±5)% for the energy carried away by the neutrinos and muons to which the FD 
is not sensitive. ^13^1^ The energy is proportional to the absolute fluorescence yield in air which 
at 293 K and 1013 hPa (337 nm band) is 5.05 ± 0.71 photons/MeV of energy depositedtlSI. The 
fluorescence yield pressure and wavelength dependency are accounted for^. By using one trig- 
gered tank in the geometry reconstruction the accuracy is improved with respect to monocular 
data (i.e data recorded only by the fluorescence detector). 



3 Energy spectrum from surface detector data 

The lateral distribution of signals is a robust measurement, the only quality criteria required for 
the surface detector data is that the station that recorded the highest signal is surrounded by 6 
active stations. This condition rejects events that might be affected by the array borders. 

The reconstructed signal at 1000 m from the shower axis on the ground level, S(IOOO), is 
a good estimator for the energy of the cosmic ray. Due to the attenuation in the atmosphere, 
S(IOOO) depends on the zenith angle: an air shower developing vertically produces a smaller 



Figure 3: (left hand side) Energy calibration. The relation between S^s and energy is almost linear and is shown 
with the continuous line, (right hand side) The fractional difference between the assigned energy and the FD 
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signal than an inclined shower produced by a cosmic ray with the same energy. The constant 
intensity method^ is exploited to obtain the zenith angle correction: it assumes that the cosmic 
ray flux is isotropic in local coordinates, i.e. the number of events above a certain threshold 
energy is constant as a function of cos^ 9. This hypothesis leads to the correction function for 
S(IOOO) shown in Fig. [2^1eft). It is a second degree polynomial S'(1000)(x) = S3s-{l + a-x + b-x) 
where x = cos^ 6 - cos^ 38°, a = 0.94 ± 0.06 and b = -1.21 ± 0.27. The quantity obtained by 
using the empirical fit shown in the same figure, S^s, represents the signal at 1000 m the very 
same shower would have produced if it had arrived from a zenith angle of 38°. This angle 
corresponds to the median of the zenith angle distribution of the SD data. The number of 
events above a certain S^s is zenith angle independent. In principle the attenuation might be 
energy dependent, because showers with higher energies develop deeper in the atmosphere and 
can be observed before their maximum development. This effect was found to be negligible. 
Distributions of the number of events for equal cos^ 6 intervals as a function of S(IOOO) and S^s 
are shown in Fig. [2] (right). After applying the correction function the zenith angle dependence 
disappears. 

The transformation from S^s to energy is obtained by so called golden hybrid events. These 
are air showers that triggered both the SD array and the FD. Only a subsample of high quality 
FD measurements is used for the energy calibrations: For example the reduced foi' the fit 
is required to be less than 2.5, data are selected only if measurements of the vertical aerosol 
optical depth are available and the maximum of the shower development is required to be in the 
field of view of the detector. Moreover, only events with a fraction of Cerenkov contribution of 
less than 50% of the total light, an uncertainty on the position of the shower maximum smaller 
than 40g/cm^ and a relative total energy uncertainty less than 20% are used in the calibration 
procedure. 

The relation between S^s and the FD energy is shown in Fig. [3] (left). As can be seen, it 
exhibits a power law correlation with a relative dispersion of 19 it 1% (Fig. [3] (right)). The 
uncertainties in the determination of both FD energy and SD signal are assigned on an event by 
event basis. The energy uncertainty 8%) includes also propagated atmosphere uncertainties 
and uncertainties from the air shower geometry reconstruction. The 5*38 uncertainty(~ 16%) 
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Figure 4: (left hand side) The differential flux as a function of energy. Also shown are the number of events in 
each bin. (right hand side) The fractional difference between the Auger vertical and HiRes 1 data and a spectrum 
with an index of 2.69^. 



contains the lateral distribution function assumption, the shower to shower fluctuation and 
reconstruction accuracies. The best fit gives the transformation from S^s to energy as E = A-S^g, 
A = 1.49±0.06(stat)±0.12(sys)[10i^eT^] and B = (1.08±0.01(stat)±(0.04)(sys)) with a reduced 
of 1.1. 

The spectrum built from 2 • 10^ events recorded until August 2007 is shown in Fig. SI 
The acceptance is computed by simple geometrical considerations and from the continuous 
monitoring of the configuration of the array The data set used for obtaining the energy 
spectrum contains only events with energies greater than 3 • 10^^ eV, since only above this energy 
the array is fully efficient. Due to reconstruction and trigger efficiency issues only events with 
a zenith angle of less than 60° are included in analysis. The integrated exposure for this period 
is 7 • 10^ km^ sr yr. The uncertainty of the acceptance is less than 5%. Having a duty cycle 
of almost 100 % the vertical spectrum from the SD has the lowest statistical and systematic 
uncertainties. 

In Fig. S] (right) the fractional difference between the vertical spectrum and a power-law 
cx is illustrated. Two spectral features are clearly visible: the so-called ankle at energies of 

~ 10^^'^ eV and a flux suppression at energies above ~ 10^^'^ eV. A continuation of the spectrum 
as a power law with index 2.69 at highest energies predicts 167 ± 3 events above 10^^'^ eV and 
35 lb 1 above 10^*^ eV, whereas we observe only 69 events and 1 event. The hypothesis of a 
pure power-law can be rejected with a significance of 6 standard deviations, independent of the 
energy scale uncertainties. 



4 The Auger spectrum: combining the hybrid and surface detector measurements 

Air showers measured by the SD array with a zenith angle between 60° and 80° are used to 
determine an independent spectrum. The procedure to derive the energy is equivalent to the 
vertical events, but instead of using 5*38 the shower size is determined from the relative distri- 
butions of the two-dimensional muon number densities at ground level ^1^. The normalization 
factor of the muon map, A'^ig, gives the total number of muons relative to a shower initiated by 
a proton with an energy of 10^^ eV. The relation between A^ig and the hybrid energy is shown 




Figure 5: (left hand side) The energy calibration for the surface detector events with zenith angle between 60° 
and 80° . (right hand side) The exposure for the three data sets used to determine the cosmic ray flux. 



in Fig. [5l The statistics is rather low compared to the vertical energy calibration, but a clear 
almost linear dependency is seen. 

The acceptance calculation for this set of events is purely geometrical and the threshold 
energy above which the trigger efficiency is more than 98% is 6.3 • 10^^ eV. Above this energy 
the integrated exposure until the end of February 2007 is 1510 km'^ sr yr; 29% of the equivalent 
acceptance for vertical events^^. 

Another data set available to measure the cosmic ray flux contains hybrid events. The 
hybrid exposure calculation relies on an accurate simulation of the fluorescence detector and 
the atmosphere. A large sample of Monte Carlo simulations are performed to reproduce the 
exact conditions of the experiment and the entire sequence of given configurations, from camera 
pixels to the combined SD-FD data taking of the observatory. The rapidly growing array, as 
well as the seasonal and instrumental effects are reproduced in the simulations within 10 min 
time intervals. The systematic uncertainty in the hybrid spectrum is currently dominated by 
the calculation of the exposure and reaches 20% in the low energy range. The advantage of 
the hybrid measurement of the energy spectrum^SH! ^j^e coverage of the energy range between 
10^^ eV and 3 • 10^^ eV. 

The exposure for the three data sets is illustrated in Fig. O The hybrid measurement 
extends to the lowest energy range. Including the use of inclined SD data improve the statistics 
noticeably as this set has about a quarter of the event totals of the vertical data set in the 
highest energy region. 

All spectra are affected by the 22% uncertainty in the FD energy scale, the main contributions 
coming from the determination of the fluorescence yield (14%). The profile reconstruction itself 
contributes with 10%. The absolute calibration of the telescopes is done every few months and 
has a contribution of 9.5%. The correction from the vertical aerosol optical depth measurements 
is 5-18%, giving an uncertainty on the energy of 4%. The advantage of the hybrid energy 
determination over the only SD assignment (e.g. as in the case of AGASA experiment) is 
that the uncertainties are experimentally driven and can be improved. A single surface detector 
experiment determines the primary cosmic ray energy with the help of simulations of air showers 
therefore the uncertainties are driven by theoretical uncertainties which are harder to decrease. 
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Figure 6: The energy spectra obtained at the Pierre Auger Observatory as presented in' — '. The combined 
spectrum is shown with filled circles. Open circles represent the energy spectrum obtained from events with a 
zenith angle of less than 60° , triangles represent the hybrid measurement and open squares the energy spectrum 
obtained from surface detector data with a zenith angle between 60° and 80° . Only statistical uncertainties are 
shown. 



The energy spectra obtained with the three methods are illustrated in Fig. [6l The systematic 
uncertainty of the energy does not affect the relative comparison of the three spectra. The 
agreement is well within the independent systematic uncertainties, the difference between the 
overall normalizations is at a level of less than 4%. 

The Auger combined differential flux, shown in the same figure multiplied with the third 
power of energy, is extending over the widest energy range possible and with minimal un- 
certainties. To deduce the spectrum a maximum likelihood method is applied taking into 
account the independent uncertainties of each measurement. The combined spectrum is 
dominated by the vertical surface detector measurement above 3EeV and by the hybrid spec- 
trum in the lower energy range. The spectral index changes from 71 = —3.30 it 0.06 to 
72 = -2.69 ± 0.02(stat) ± 0.06(sys) at log{E^^i,iJeY) = 18.65 ± 0.04, and above lO^^-^eV 
to 73 = -4.2±0.4(stat). 

5 Conclusions 

The energy spectrum has been measured at the Pierre Auger Observatory with three independent 
data sets and the agreement is better than 4%. 

A flux suppression at the highest energies has been established with a significance of 6 
standard deviations. Combined with the recent result that the highest energy cosmic rays are 
anisotropic and a good correlation only with relatively nearby sources has been found , it hints 
towards a GZK effect. 

A change of the spectral index occurs at about 4 EeV which might indicate either the tran- 
sition from galactic to extragalactic origin of cosmic rays or a propagation effect. To distinguish 

93' 

between the models a determination of the mass of the cosmic rays over the entire energy 
range is necessary. The exact location of the spectral index change and an accurate measure- 
ment of the spectrum below 3 EeV will be possible in the near future with AMIGA and HEAT 
extensions of the observatory. 

The exact shape of the energy spectrum in the highest energy range, which can be used to 
distinguish between acceleration scenarios'^U, will be determined with greatly improved precision 
over the next 10 years of data taking. 
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